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1. INTRODUCTION 

While the existence of Dark Matter (DM) has now been well established, it is still unknown 
what object, or objects, makes up this dominant component of the matter in the universe. Perhaps 
the best motivated candidate is a new neutral particle whose mass is around the weak scale (0(1 — 
1000) GeV) and whose couplings to Standard Model (SM) fields are somewhat below the weak 
scale. So far direct detection experiments have seen no clear signal consistent with DM recoiling 
against atoms [IH3] , with several notable exceptions [IHS] • At the same time the LHC is performing 
a multitude of searches, many of which require substantial missing energy. Although these searches 
are sensitive to the production of DM, most are designed to search for models where the DM is 
made in the cascade decay of a colored particle, e.g. supersymmetric searches, and are not generally 
applicable to all models. 

Recently, it has been noted that there is a class of more model independent searches for DM 
production that can be carried out at the LHC 1 , based around the production of DM in association 
with monojets [9rll5|. monophotons |13|. I15j . mono-VF [IB] and also a more inclusive multijet 
search [T7j. This class of searches can place strong constraints on the properties of DM that are 
complementary to those from direct detection searches. In particular, the LHC does not suffer 
from a low mass threshold, nor is the spin-dependent (SD) bound considerably weaker than the 
spin-independent (SI). Although the collider SI bounds for DM mass above ~ 10 GeV are weaker 
than direct detection, the SD bounds are often stronger. In addition, the DM is being produced by 
the experiment rather than relying on a galactic component in our vicinity and so is not sensitive 
to unknown astrophysics. 

A dedicated DM search, using the shape of the leading jet transverse momentum, px, distribu- 
tion has been carried out at CDF on 6.7 fb _1 of data |18j . CMS and ATLAS have both carried out 
cut-and-count based monojet [H2 [20] and monophoton [2H [22] analyses using ~ 5 fb _1 of data. 
Given the copious rate for the production of Z's in association with jets the current round of exper- 
imental results are already dominated by systematic (rather than statistical) forms of uncertainty. 
In the near future this will also be the case for the smaller (but by no means small) production rate 
of Z in association with a photon. Given that the analyses are therefore systematics limited the 
theoretical community should investigate the possibilities in which the systematic errors associated 
with the theoretical predictions can be reduced. 

The most obvious mechanism to reduce the theoretical systematic uncertainty is to provide Next- 
to-Leading Order (NLO) predictions for the irreducible (and reducible) backgrounds. However, 
NLO predictions for Z + (j/j) have been available for over a decade [23} I24j . In addition, recent 
developments in matching showers to NLO predictions have resulted in a publicly released matched 
shower prediction for the Zj process within the POWHEG-BOX formalism [25]. Therefore, until 
the completion of the NNLO Z + (j/'j) cross sections potential improvements in the theoretical 
predictions for the dominant backgrounds are limited (although very recently, EW corrections to 
monojet production have been computed [26J). 

The remaining scope for reducing the theoretical systematic errors resides in improvements in 
the modeling of the signal process. Thus far experimental analyses have relied on Leading Order 
(LO) predictions for the signal and as a result are exposed to the inherent issues associated with 
a LO prediction, namely a large uncertainty in rate and shape (the shape issue can be improved if 
matched shower predictions are used). By improving the theoretical predictions for DM production 
by including NLO corrections one therefore reduces the overall scale dependence (reducing the rate 
uncertainty) and typically obtains a larger value for the cross section. Therefore, a limit obtained 



For constraints on DM coming from other colliders see [3 [8]. 
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using a NLO prediction for the total rate is both stronger (since the cross sections are larger) and 
more accurate (since the rate uncertainty is reduced). 

Furthermore, the typical monojet search cuts on the missing transverse momentum, pj>, and 
the leading jet pt are usually mismatched and the searches employed at the LHC are more inclusive 
than the monojet names suggests, allowing up to two jets in the event. Thus, there is a region 
of phase space that has two jets whose individual pr is below the $ T cut, but above the jet pt 
requirement, which together recoil against the DM, that is not included at LO. This contribution is 
naturally included in predictions which have access to the DM plus two-jet phase space, e.g. NLO 
and matched shower predictions. 

For these reasons we therefore implement the production of a pair (of fermionic) DM particles in 
association with either at jet or a photon into the MCFM |27H29j code, which is available publicly. 
We will study a range of phenomenologically interesting operators, primarily (but not restricted 
to) the effective theory in which the particle responsible for mediating the DM/SM interaction is 
very heavy. We will present our results in a format which allows easy recycling of our results, 
more specifically we present the results for the SM and DM parts of the calculation in a factorized 
manner. This allows easy implementation of our results to other new physics scenarios which 
produce monojet or monophoton signatures. 

This paper proceeds as follows, in section [2] we provide the necessary DM and SM ingredients 
to construct our NLO calculation. In section [3] we investigate the NLO phenomenology associated 
with the monojet final state. Section [4] presents a similar study for the monophoton channel. We 
draw our conclusions in section[5} In appendix|A]we include some spinor definitions and appendix [B] 
catalogues some of the formulae obtained during our calculation. 



The strategy we will follow to provide NLO predictions for monophoton and monojet DM 
processes is to factorize the problem into SM production and followed by a BSM "decay" in the 
final state. That is, we will consider operators of the form OsmOdm where the particle coupling 
the SM to the dark sector is exchanged in the s-channel. For instance, the amplitude for monojet 
production through the vector operator is of the form, 



Situations involving t-channel mediators, for instance squark exchange, can still be considered 
by carrying out a Fierz transformation. These will then involve contributions from multiple s- 
channel operators. Our implementation in MCFM allows for this full generality but we do not 
consider such a case here. We will be predominantly interested in DM production which proceeds 
via an effective field theory (EFT), where the exchanged particle is integrated out. We we also 
briefly provide examples for the case of a light mediator for a subset of our operators. The effective 
theory is well motivated provided that the mediating particles are heavy (> few TeV), however 
whether or not the mediating propagator is included bears little impact on the NLO calculations, 
our MCFM implementation can calculate either in the EFT or the full theory. The operators 
we consider that involve SM quarks correspond to vector, axial-vector, scalar and pseudo-scalar 
exchange, and we consider one operator that couples gluons to DM, 
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NLO CALCULATIONS OF DM PRODUCTION 




(1) 



(X7MX)(g7 M g) 
A 2 

(X7^75X)(g7 M 75g) 
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n _ (xx){Gg U G a ^ l/ ) 

u 9 — a s ^3 ; W 

O s = p , (5) 

P5 = A3 ' ( ' 

We restrict our focus to the above operators, which provide a representative sample of the phe- 
nomenologically interesting models. Due to our factorization approach, more general cases, say 
for instance an operator of the form (xifj,x)(Ql^l5 ( l): can De readily obtained from the results 
we will present. The operator Oy(A) has a simple UV completion involving exchange of a vector 
(axial-vector) boson, of mass M and width V. The full theory corresponds to the replacement 
A 2 — > (s^x — M 2 + iMY) / g x g q . Operator O g is induced at the loop level and the simplest UV com- 
pletion involves a heavy scalar and heavy fermions. The scale of the new physics, M, is typically 
lower than for 0yM) an d so the effective theory has a smaller range of validity. 

For the scalar and pseudo-scalar operators we have written the couplings as scaling with quark 
mass. This is what is expected if an assumption of minimal flavour violation (MFV) is made, in 
which case the only flavour violating spurions are the Yukawa matrices. With SU{2) invariance 
requiring an implicit Higgs field insertion the operators then scale with quark mass. Here we only 
consider the flavour diagonal part of these operators, flavour violation in DM couplings leads to 
other interesting signals |3Q|, 131] . If the MFV assumption is loosened one has to contend with 
strong constraints from flavour observables. One possibility would be if the DM couplings were 
0(1) in the quark mass eigenstate basis, although such a model would be highly tuned. Our 
MCFM implementation is sufficiently flexible to allow either possibility. We focus in this paper on 
the more motivated MFV case and present results for that. Because of PDFs and the suppressed 
light quark couplings, there is a difference between the behaviour for the case of DM coupling to 
the first five generations of quarks and to the top quark [32J. Therefore we will address the light 
quarks (u,...,b) and top separately. 

Finally, unlike in the effective theory, i-channel operators present an additional problem upon 
UV completion. Since the exchanged state has to be coloured there are additional diagrams that 
must be included in order to achieve NLO accuracy. Therefore, for t-channel operators our NLO 
results can only be used for the case where the mediator has been integrated out. 

In order to facilitate the calculation we will simplify the problem further by calculating helicity 
amplitudes. For the massless SM production amplitudes this introduces a dramatic simplification 
in the number of independent calculations, since helicity amplitudes naturally involve projections 
of the form (1 ±7 5 ). This allows us to determine the results for the various operators from common 
building blocks. Since the DM particles are massive care must be taken because the mass spoils 
the chiral symmetry. 



2.1. Dark Currents 

We wish to use helicity methods in order to calculate our dark currents. However, since the 
fermionic DM is massive helicity is not a good quantum number. Nevertheless, massive fermions 
must satisfy the usual sum rule when summing over polarization states, 



^Usip, m)u s (p,m) = f + m (7) 
v s (p, m)v s (p,m) = f-m. (8) 



s=± 
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Therefore, one can re- write a massive spinor in terms of two massless spinors [33] provided that 
the completeness relation above is preserved. In order to write massive spinor bilinears in terms of 
massless ones we use the prescription of refs. \3A\ 155] rewriting momenta (p\ and p2 with pf = m 2 ) 
for two massive particles, of mass m, in terms of two massless particles (k\ and /c2, with kf = 0) 
such that, 

Pi = X A^K + (9) 
= + l -^K (io) 

where, j3 = yl — 4m 2 /si2, and su = (pi +P2) 2 = 2k\ • /c2- This decomposition has the advantage 
that p\ + P2 = ki + k2- The helicity states for the massive spinors (u, v, etc) are defined in terms 
of the massless spinors |1+), |2+), through the following relations, 

£-1/2 ^-1/2 
u±(p2,m)= ^ 2± j {l T \ty 2 + m), v±(p 1 ,m) = T^±j if 1 ~ m)\2 ) , (11) 

where we have introduced variables (3± = ^(1 ± (3). A list of spinor definitions is provided in 
appendix [Aj and we encourage the interested reader to inspect ref. |35j for additional properties 
and relations of massive spinors. We are now in a position to define the currents we will need for 
our study, namely those involving the insertion of a 7 matrix, 

V(2± if) = u±( P 2)rv T ( Pl ) = (2 ± | 7 H1 T ) (12) 
V^(2+, 1+) = V+faWv+fa) = 2—(h - k 2 r , (13) 

777 

V(2- I") = «_(p 2 ) 7 ^_(pi) = 2^(fc 1 -A !2 )''. (14) 

Note that the helicity conserving currents are identical to their massless counterparts, and that the 
helicity violating currents vanish in the massless limit as required. For the cases in which the dark 
matter is mediated through the exchange of an axial vector we will need the following currents, 

V£(2± if) = ^feh^Mpi) = T(l - 2f3 + )(2 ± \r\l T ) , (15) 



Tft 

V 5 "(2+ 1+) = u + (p 2 h^v + ( Pl ) = -2^(h + k 2 r , (16) 

711 

V£(2~ 1") = u„( P 2h^v-( Pl ) = +2pj(A;i + k 2 Y • (17) 

In addition to the currents listed above we will also need the following scalar currents, 

5(2+ If) = u ± (p2)v T { Pl ) = , (18) 
5(2+ 1+) = u + (p 2 )v + ( Pl ) = (2/3+ - 1) [21] , (19) 
5(2", I") = u-( P 2)v-( Pl ) = (2/3+ - 1) (21) . (20) 

Finally we will consider the decays of dark matter through a pseudo-scalar current 

5 5 (2± If) = u±(p2h 5 v T ( Pl ) = , (21) 
5 5 (2+ 1+) = u + (p2h 5 v + ( Pl ) = - [21] , (22) 
5 5 (2-,l=) = u_(p2)7st;_(pi) = (21) . (23) 

Note that in all cases, the m — > limit is clearly reproduced correctly. 
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2.2. SM Vector Currents 

We now define the currents needed to describe the SM production of a gluon and a vector current, 
Oy. Such results are already present in the literature [231 ES]> namely amplitudes involving the 
production of a Z (or more precisely a virtual photon 7*) and jets. In order to efficiently recycle 
these amplitudes one must first remove the unwanted decay of the Z into two massless leptons. 
The general strategy is thus to re-write the amplitude for the Z in the form, 

2+ 3r,4+ 5r) = 4°^(l+ 2+ 3-) x J*(4+ 5r). (24) 

The current can then be replaced by the DM current of choice to produce the desired result. 
Similar techniques were used in ref. [37] for the calculation of the amplitudes for Z^/j in which the 
photon was radiated from the final state leptons. 

As an example we present the amplitudes for the tree-level production of a jet (or photon) in 
association with a vector operator, with momentum P flowing in the XX system, 

2 + 3^ _ 1 (3|(1 + 2)|^|3) 

4°^(l+2-,3r) = ^(1+ 2+ 3T) ((ab) o [ab] , (1 o 3)) . (26) 

Since the extraction of the currents needed for the NLO corrections are now determined in terms of 
the known results in the literature |23|. 136] we refrain from writing them explicitly here. However, 
since they are useful building blocks for other mono jet searches we present a full list of the currents 
(for virtual and real corrections) in Appendix [Bj Before moving on to discuss other SM currents we 
note that the monophoton amplitudes are naturally related to the pieces of the monojet amplitudes 
which are subleading in colour. 



2.3. SM Scalar Currents 

In this section we provide the helicity amplitudes for the "SM" production of a scalar in asso- 
ciation with a jet or a photon. In the previous section we described the extraction of the vector 
currents from the existing literature results for Z/^f*+ jet. However this calculation has no obvious 
analog in the SM since the Higgs couples to massive fermions and we consider five massless flavours. 
Therefore, we perform the calculation directly. The NLO calculation of DM plus monophoton for 
the scalar operator has been studied previously, using traditional matrix element techniques, in 
|38| . We will focus on the case of scalar and pseudo-scalar couplings which are proportional to 
mass, similar to the SM Higgs couplings and therefore we could have used the H + b results from 
the literature. However we wish to remain general enough to allow for couplings independent of 
mass, therefore we drop the mass terms in the SM production which are not part of the Yukawa 
coupling. Since (to the best of our knowledge) the helicity amplitudes for these processes have not 
been written down before, we present them in full in this section. 

For completeness we note that the case of scalar coupling to the top quark is particularly 
interesting since there is no tree-level monojet diagram and the LO result is due to a loop of top 
quarks [32], these results can be easily extracted from the Higgs plus jet results in MCFM [39]. 

We begin by listing the two independent tree-level amplitudes, 

4°>(w,3 3 + ) = <|^, <m 

^C"'W = 5Pi' (28) 
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The virtual corrections to this amplitude consist of leading and sub-leading colour contributions. 
The analytic forms of these expressions are rather simple, primarily because the insertion of the 
scalar operator does not increase the tensor rank of any loop diagrams. We calculate the four- 
dimensional cut-constructible pieces using the quadruple cut technique [10]. We have checked 
our virtual results against a numerical implementation of D-dimensional unitarity [H], finding 
agreement. This check also confirms the lack of bubble pieces from the expressions. The expressions 
for the virtual amplitudes are very simple, the leading colour pieces have the following form, 

J8»K,KX) = orAf(-U(^-) e +(^-)')^J^,=^-)) , (29) 



V 6 \\— s 12j \—S23/ J \— S123 — S123 

. Sl3 + «23 / Qn N 

+Cr 2[23]M • (30) 

The subleading in colour amplitudes have the following form (these are also the monophoton 
amplitudes) 

+ U_J=^,=^-)\ , (31) 

\-Sl23 -5123/ / 



|T / -523 -S12\\ . 513 + 523 
\-5l23 -5123// 2 [32J [13J 

Unlike the vector case of Appendix [Bl there is no Ward identity for the scalar operator. This results 
in a further UV counter term in addition to those specified in Appendix [Bj For example, in the 



MS-scheme one must include the counterterm — cpJ-.A;? in (29)-(32). For the calculation of the 



real corrections we will need the following amplitudes which involve the emission of an additional 
parton with respect to the Born monojet process 

4 0) (l-,2r, 3 -,4 3 ) = [14] . (33) 

Af(l-,2=,3+, 4+) = {u) { {2 } 3) {u) , (34) 



4 (0) n - 2 - 3- _ <23>[42]<1|(2 + 3)|4] (13) 2 (2|(1 + 3)|4] 
S Kq ' ~ q ' 9 ' 9 ] ' 5 2 34 (34) [23] [34] + s 1S4 (14) (34) [34] 

(13)(1|(2 + 3)|4] 



(35) 



(14) (34) [23] [34] ' 

4 0) (l-,2r,3+ 4-) = 4 0) (l-,2-,4-,3+) . (36) 

In addition we will need the four-quark amplitude which is given by, 

,(0) n - 9 - „_ _ <13)(2|(1 + 3)|4] (23)(1|(2 + 3)|4] 
A 8 V 9 ,2,,d Q ,fy- si34 (34)[43] + S234 ( 34 )[43] • {S7) 

Finally for the case involving monophoton production the amplitudes we need are as follows, 

A (°)n- 9- o- a-\ . g 1234 [12] 

A S U, >*q,*g,*y) ~ [u] [23] [13] [24] > 
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^(^g ' j 3 ff , 4^ ) 



(12) 



(14) (13) (23) (24) 
<14> <2|(1 + 4)|3] 
S134 (13) [41] 

(12) S124 



(24)(1|(2 + 4)|3] 
S234 (23) [42] 



+ 



A 



(13) (23) [41] [42] 
(0) f „ 



4")(3o4) 



(39) 

(40) 
(41) 



2.4. SM Gluonic currents 



The production of DM in association with a jet through the gluon operator 

O g = a s m{G ^ Ga ^ ] , (42) 

is simple to implement due the large existing literature associated with Higgs boson production 
in the heavy-top effective theory. In fact, modulo changes in the coupling the process is identical 
to H — > (66)+jet. Therefore we simply modify the existing MCFM routines accordingly [39, 42J. 
Clearly no monophoton signature is produced from this operator. 



3. MONOJET PHENOMENOLOGY 

In this section we present some phenomenological studies using the results derived in the previ- 
ous sections. We have implemented our NLO calculations into MCFM, taking full advantage of the 
codes existing architecture to obtain our predictions. More specifically, we use the dipole subtrac- 
tion scheme of Catani and Seymour [33] to render both the virtual and real corrections separately 
finite. We use the following MCFM default electroweak (EW) parameters in our calculation, 

M z = 91.1876 GeV , M w = 80.398 GeV , 
T z = 2.4952 GeV , T w = 2.1054 GeV , 
G F = 0.116639 x 1(T 4 GeV~ 2 , m t = 173.2 GeV . 

The remaining EW parameters are defined using the above as input parameters. In order to avoid 
DM being charged under SU{2) the scalar operator proceeds through a Yukawa coupling and is 
thus proportional to the quark masses. In these cases we use m c = 1.5, nib = 4.7 GeV and m t as 
defined above with all other quarks kept massless. For the light quarks (m c and nib) the masses 
are only retained in the form of a Yukawa coupling, and are not retained in the kinematics of the 
matrix element. For a pure monojet analysis this is of little consequence. In the case of the top 
quark, which can contribute to monojet processes at the loop level, the mass is fully retained. Our 
default PDF choice is CTEQ6L1 for LO and CT10 for NLO calculations gl]. 

In this section we will use the effective theory prescription defined earlier. These EFTs are 
adequate descriptions of a more UV complete model, provided the mediator scale A is large enough. 
Our implementation of the various operators into MCFM allows for flavour dependent couplings 
to be used, so that up- and down-type operators can be studied individually if so required. In 
these examples we will consider the operators to be flavour diagonal, coupling to each flavour with 
the same coupling, taken to be 1. The only exception to this rule are the operators Os and Ops, 
which as discussed above couple with quark-mass dependent couplings. 
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3.1. Monojet inclusive cross sections 

We now turn to phenomenology, we begin by considering the production of a pair of DM particles 
in association with a jet. In this section we focus on the LHC operating at 7 TeV and we will 
use cuts inspired by the recent ATLAS and CMS publications [TD1 [2D] requiring events to pass the 
following phase space cuts, 

$ T > 350 GeV , j4- > 100 GeV ' < 2 > ^h,h < 2 - 5 • ( 43 ) 

Jets are defined using the anti-fey algorithm with R = 0.4. Experimental searches also implement 
a veto on more than two jets, however this kinematic configuration is not covered by the NLO 
calculation and as such this cut is meaningless in our studies. As discussed earlier, at LO the 
mismatch in the $x and jet-py requirements means that only events with the jet pt above the j£y 
cut will be included. At NLO there is a new region of phase space with a jet and an additional 
parton (which may or may not pass the jet algorithm) each with a px below the $t cu t> that can 
conspire to produce enough $ T so that the event passes all cuts. We note that the cut on the 
azimuthal separation between the two jets A4>j 1 j 2 only enters our calculation at NLO and as such 
is a fairly weak cut. 

Using these cuts we begin by studying total inclusive cross sections and their inherent scale 
dependence at LO and NLO. Our results are summarized in Fig. [TJ We have chosen the values 
of A for each operator to get close to the rate the experiments are presently placing a bound on, 
O(10 — 100) fb. We have made no attempt to actually place limits directly on A. In addition we 
provide the ratio of NLO to LO cross sections (the A-factor) for each operator in Fig. [2j The 
dependence on the renormalization and factorization scales, taken to be equal (f/,R = fj,p = fj,), is 
illustrated by the shaded band linking the predictions obtained at fj, = 2m xx and \i = l/2m xx , 
whilst the central scale choice \i = m xx is illustrated by the curve inside the shaded band. From the 
various curves it is clear that the NLO corrections to the total cross section are sizeable, and vary 
operator to operator. As is natural by going to NLO the overall scale dependence is smaller than 
the LO prediction. Before we comment on the results for each DM operator individually we observe 
that all the plots have a similar shape. Namely a flat plateau in the light m x region in which the 
cross section is insensitive to m x , and then a rapid fall off at a DM mass of around 100 GeV. These 
features have been observed in previous phenomenological studies [9TH5| and indeed are clear in 
the experimental constraints on A \19\ 120] . For light DM the analysis cuts are sufficiently hard 
that a \/s = 7 TeV collider sees no difference in terms of phase space restrictions in producing DM 
with mass 10 GeV or 0.1 GeV. This explains the plateau in production cross section as a function 
of DM mass. Then as the DM mass becomes a non-negligible scale in the process the phase space 
restrictions begin to rapidly drive the cross section downwards, resulting in the distinctive fall off 
in each of the plots. The point at which this fall off begins depends on the exact mass dependence 
of the cross section and is therefore operator specific. As the operating energy increases, the fall off 
moves to larger m x , for every operator. As such, the limits on A obtained using the 8 TeV data set 
should show a considerable improvement (beyond the simple y/s rescaling) for DM masses around 
0.1-1 TeV. 

We now examine each operator in detail, firstly we observe that the vector and axial-vector DM 
operators show similar behaviour in terms of AT- factors and scale dependence. This is unsurprising 
since in the massless limit the only terms which are sensitive to the axial nature of the coupling 
are the four-quark amplitudes, which are a small part of the total NLO cross section. As a result 
it is clear that as m x — > the results must be similar. As the the DM mass grows, so does the 
difference between the operators, with the axial operator being smaller than the vector over the 
entire mass range. The scale variation for these operators is also similar, with typical LO values 
around ±20% and NLO values around ±10%. The A-factor for both these operators is around 
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Figure 1: LO and NLO cross sections for DM production in association with a jet at the 7 TeV LHC. The 
solid line indicates the cross section obtained with the default scale (i = m XXl the shaded band represents 
the deviation from this scale when the scales are varied by a factor of two in each direction. The phase 
space cuts described in the text (43) have been applied. 
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Figure 2: if-factors for DM production in association with a jet at the 7 TcV LHC. The solid line indicates 



the i-L-factor obtained with the default scale fi 



the shaded band represents the Jf-factors obtained 



using scales varied by a factor of two in each direction. The phase space cuts described in the text (43) have 
been applied. 



1.4. Since, for this operator, the cross section scales as A -4 one naively expects an improvement 
on the limits on A of around a ~ 10% if the NLO rate were used compared to the LO one. 

The O g operator shows some striking differences with respect to the Oy and Oa operators. 
Firstly the iv-factors are much larger, with values between 2 and 2.5, indicating much larger NLO 
corrections for this operator. Secondly, the scale dependence is large at both orders. The LO scale 
dependence is around 50% reducing to around 25% at NLO. These effects are reminiscent of the 
Higgs effective operator, in which large scale dependence and X-factors are common. However, 
there is a crucial difference between the two calculations which should not be overlooked, for 
the case of a light Higgs the final state phase space is extremely restricted by the narrow s- 
channel resonance. In our effective theory, no such restriction applies, allowing for the possibility 
of enhancing the NLO effects by sampling over a larger phase space. Since the operator scales as 
A~ 3 using a NLO cross section to set a limit on A should result in improvements of the order of 
~ 16%. 

The scalar and pseudo-scalar operators (for m& and m c couplings) (0^'^ and Opp) behave 
in a rather different manner to the other operators. They are naturally suppressed by the factor 
m q /A compared to the vector, and axial cases, which results in much smaller cross sections. For 
this reason we have chosen to use a smaller value of A in our analysis of these operators (50 GeV 
compared to 500 GeV for the other cases). The LO diagrams require the presence of a heavy sea 
quark (charm or bottom) in order to be non-zero, therefore it is natural to expect a large i^-factor 
at NLO arising from contributions in which the process is initiated by two gluons. This is indeed 
the case, we observe a /f-factor of around 2-2.5, leading to an improvement in the A bound of 
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Figure 3: NLO predictions for the missing transverse momentum spectrum for signal and background 
obtained using MCFM. We present results for two operators (Oa and O g .) The DM mass is 100 GeV, and 
the scales have been chosen to be equal to the DM invariant mass [i — to xx . 



~ 25%. The scale dependence is markedly reduced at NLO compared to LO, (around 5% at NLO 
compared to around 20% at LO). However, given the large K- factor going from LO to NLO it is 
unlikely that scale variation alone gives a reliable estimation of theoretical uncertainty (even at 
NLO) for this operator. 

Finally, we consider the production of a monojet signature using the Og operator. The en- 
hancement which comes from including the loop induced processes in which the DM couples directly 
to the top quarks, is very large [32] ■ We show the LO cross sections obtained using MCFM where 
for this operator we have set A = 150. The cross sections obtained using the top-induced operator 
is around two orders of magnitude larger than those for the light quarks. As a result if this opera- 
tor was used in place of the light quark operator the corresponding limit on would increase A by 
between 200 and 300% [32] . Since this processes is LO (even though it is a loop induced process) 
the scale dependence is very large, around (40-50%). This means that a sensible strategy experi- 
mentally may be to set limits on the operators induced by top and bottom couplings individually. 
Although the bottom induced A will be much smaller than that for the top induced coupling it will 
suffer from much smaller systematic uncertainties on the theoretical side, since the NLO prediction 
can be used. However, we stress that the resulting A limit obtained by either operator is unlikely 
to be within the realm of validity for the effective theory. For the top-scalar operator the limit is 
salvageable by switching to the full theory as we will investigate in the next section. The rate for 
the light-quark operators are suppressed by quark mass and so the EFT approach is brought into 
question. In a simple UV completion of a scalar mixing with the Higgs, discussed below, the rate is 
also expected to be very small. However if a, naively, tuned model is considered in which the DM 
couplings are diagonal in the quark mass basis the couplings could instead all be 0(1) and these 
operators would have non-trivial bounds. Given that the production of DM through the scalar 
operators probes a different set of PDFs from the vector operators it is still worth investigating 
the scalar operator in general. 



3.2. Monojet differential distributions 



So far we have focussed our attention on predictions for total inclusive cross sections at NLO. It is 
also interesting to consider the NLO corrections to important distributions used in the experimental 
searches. As an example of such a distribution we consider the missing transverse momentum 
spectrum For simplicity we focus on two operators (Oa and Og) for a fixed DM mass of 
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Figure 4: The differential if -factor, (daNLo/dtr)/ ' {d<jLo/d$ T ), for monojet production via either an axial 
or gluon-induced operator for m x = 100 GeV. 



100 GeV. We also use MCFM to obtain the spectrum for the dominant background contribution 
Z — > (vv)j, our results are shown in Fig. |3| As was the case at LO the $i T spectrum for the 
DM signal is noticeably harder than for the that of the dominant background Zj. This can be 
understood by the scaling properties of the two spectrums, the signal to background ratio scales 
as A 4 /sp^. Therefore, in the limit of large MET (large Sp u ), the EFT produces a harder spectrum 
than the background. This effect is dominated by the differences between the full and effective 
theory and as such is unchanged at NLO. In order to quantify the differences between LO and 
NLO for the signal distribution we present the differential if-factor in in Fig. |4| The choice 
of a dynamic scale at NLO results in a fairly stable differential if-factor with only a small growth 
with increasing 



3.3. Light Mediators: Full theory considerations. 

The results presented above were obtained in the effective theory in which the particle re- 
sponsible for mediating the interaction between the SM and the DM was taken to be heavy and 
was integrated out. It is interesting to consider the accuracy, the phenomenological validity, and 
the implications if this assumption breaks down. In particular, we note that, for the scalar and 
pseudo-scalar operators, in order to achieve DM production cross section of the size necessary to 
be observed at the LHC we must consider A of 0(50 — 200) GeV, see Fig{T] Given the typical scales 
involved at the LHC, A of this order cannot correspond to a mediator that can reasonably taken 
to be heavy, without entering the realm of strong coupling. As a result, we investigate the monojet 
phenomenology by reinstating the s-channel mediator. For simplicity we focus on two simple UV 
completions, the first being in which the mediator is a massive gauge boson which couples axially to 
the SM and DM fermions (i.e. a straightforward UV completion of Qa). The scale in the effective 
operator Q is related to the scale of the mediator by A = M§ j ^/g q g x where the coupling to quarks 
(DM) is g q (g x ). 

Secondly, we consider a UV completion of the scalar operator for top quarks only, xxtt- A 
simple example of such a UV completion is to consider a new heavy singlet scalar $ that mixes 
with the Higgs through a fiH 2 <$> operator. This leads, after electroweak symmetry breaking, to a 
coupling of the top to the mostly <3? mass eigenstate that is proportional to the top Yukawa. If this 
heavy singlet also couples to DM as g x &XX then the amplitude for DM coupling rescales as 

_> Ml (44 ) 

A3 ^ m\ x - Ml + iM$F$ ■ 1 ' 

Thus, the relationship between the full theory and the effective theory is more complicated and 
A = (M V^ X )V3. 

Results for these two theories are show in Fig. [5j we have used the same final state phase space 
cuts as in the previous section and to aid in comparisons to the EFT we have fixed g x = gt = g q = 1 , 
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Figure 5: Comparison between the full and effective field theories for axial and top-scalar operators, for DM 
mass of 100 GeV. The straight line indicates the results obtained in the EFT, shaded region indicates full 
theory results as a function of the width between the values (lower) and M$/(&tt) (upper). We have 

used the same phase space cuts as in the previous sections. The scale choice is [i = m xx . 

relating the scales of the operators in the EFT approach to the scale of the mediator as above. 
We then proceed to calculate cross sections as a function of in the full and effective theories. 
For the full theory one must also specify the width of the mediating particle T^. We choose two 
example widths, r<j> = M^/3 and T$ = M^f^ir). The results are naturally dependent on the 
choice of DM mass, for simplicity we have chosen to focus on a single value m x = 100 GeV. 

The two plots for the two different operators shown in Fig. [5] have some generic features which 
we explain first before mentioning some operator specific phenomenology. Firstly, it is clear that 
full theory asymptotes to the effective theory at large as required. On our log-log plot the 
EFT results then possess simple scaling as a function of A, and therefore (in the axial case) 

2/3 

and (in the scalar case). All of the features of the FT are dominated by the inclusion of the 

propagator in the cross section, and in particular whether or not the propagator is able to provide 
resonant enhancement to the cross section. If the propagator is able to go on-shell then the cross 
section is enhanced beyond the EFT approximation, whereas if the propagator is forced into the 
off-shell region the cross section is suppressed (dramatically for light mediators) relative to the 
EFT. 

It is simple to determine whether the resonant enhancement will be included as a function of 
m x and $f^ m , since, in order to achieve the on-shell condition one must have, 

s xx ~ Ml =► {$T n ? + Am 2 x < Ml . (45) 

Therefore if the cut is too hard, making p T too large, or m x is too heavy then the mediator 
cannot go on-shell, suppressing the cross section. For our setup this occurs at roughly 400 GeV. 
This simple kinematic argument however, has big implications for experimental searches in the 
full theory. One should endeavour to not cut away the region in which the signal peaks, therefore 
adjusting the minimum $j> cut to become a function of M^ should provide a natural way to 
optimize signal over background. 

We now make some operator specific statements, firstly we note that the limits obtained on 
A for the axial coupling are of the order 0(1) TeV, and therefore vindicate the use of the EFT 
in these searches, however the results on A are also bounded by below, since the EFT cannot in 
general be trusted for A < 200 GeV. Secondly we note that the limits on A from the top-scalar 
operator are of the order A ~ 150 GeV [32] . This is exactly in the region of M^ ,g,A phase space in 
which the EFT begins to breakdown, however, if the FT was used with slightly softer $x cuts one 



15 



1 


50 


1 


45 






11 


40 


_l 






35 


1 


30 


1 


25 



r>=M„/3 
m K =100 GeV 



10 50 100 500 1000 

MAGeV) 



5000 



Figure 6: The ratio of NLO to LO cross sections for the axial operator in the FT, for DM mass of 100 GeV. 
The red curve indicates the results obtained with width set to /3 and the blue curve indicates the results 
obtained with a width of M t f,/(8n). 



obtains cross sections in which the EFT and FT are similar. It is clear that when setting limits on 
the top-scalar operator, the FT effects should be investigated, particularly since in this region the 
validity of the EFT is correlated with the exact phase space cuts applied. Note, however, that over 
a large range of mediator mass the EFT actually underestimates the size of the DM production 
cross section, leading to the present LHC constraints, carried out in the EFT, being conservative. 
We note that the limits set on A by the charm- and bottom-scalar operators are in the region where 
the EFT is extremely dubious, since the cross section here is suppressed by around five orders of 
magnitude. 

Although thus far we have only considered fixed widths, r,*, for very wide mediators a running 
width, r(s), may be more appropriate [35]. This corresponds to the replacement in the mediator 
propagator of, 

(*xx " M l) + iM^t -> (sxx ~ M l) + *v^c f {s xx ) ■ (46) 

When the mediator is heavier than all the particles it couples to, the running width may be 
approximated by, T(s) = ^/sT^/M^. Given the flexibility of our code it is straightforward to 
implement the above changes to include running width effects. As an example, for the top scalar 
coupling and a mediator mass of 1 TeV we find the difference in the total cross section between 
the fixed and running width is sub-percent for the narrower case. For the wide mediator, the effect 
is larger, as expected. The running width increases the cross section by ~ 15%. 

Finally, it is interesting to consider the impact of our NLO corrections on the axial operator in 
the full theory. Therefore we plot the ratio of NLO to LO cross sections as a function of for 
our two width choices in Fig. [6j There is a clear dependence on M^, with the iT-factor growing as 
a function of Ma,. This dependence is controlled by the width of the mediator and the resulting 
shape of the Breit-Wigner. As expected for heavy mediators the effective theory inclusive .fT-factor 
is obtained. 



4. MONOPHOTON 

In this section we describe the phenomenology associated with DM production in association 
with a monophoton. In order to define photons in a hadronic environment an isolation criteria must 
be applied. This isolation reduces the contributions from the production of secondary photons in 
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the decays of certain types of hadrons. Typically an experimental isolation criterion requires the 
amount of hadronic energy inside a cone around the photon to be less than a fixed input, i.e. 

E^ ad <E^ mx with R = ^Acf)" 2 + At? 2 . (47) 

had£i?o 

We use typical values at the LHC for the isolation cuts, i.e. Ro = 0.5 and E^ 8 * = 5 GeV. On 
the theoretical side this form of isolation introduces complications at NLO since photons radiated, 
through bremsstrahlung, from final state fermions induce a collinear singularity. This singular- 
ity has no corresponding singularity on the virtual side. In order to render the calculation finite 
the contributions from the fragmentation functions must be included. These fragmentation func- 
tions require non-perturbative input needed to provide boundary information at a certain scale. 
Our treatment of the fragmentation contributions is the same as in previous photon processes 
implemented in MCFM [29, 37J, and we use the fragmentation functions of [46J. The vector and 
axial- vector production of DM along with a photon have been recently considered [17], although 
the effects of photon fragmentation were not included. 

4.1. Phenomenology 

The procedure we follow for the monophoton examples we consider is very similar to that 
described earlier, section [3j We base our cuts around those used in current experimental analy- 
ses [2U [22] , more specifically we require our final state to satisfy, 

$ T > 140 GeV , pi > 150 GeV , |?? 7 | < 2 . (48) 

We present LO and NLO cross sections, under these cuts for various operators in Fig. [JJ K- 
factors obtained for these operators are illustrated in Fig. [8} We observe that, as was the case for 
monojet examples, the axial and vector results are similar to each other, and the scalar and pseudo- 
scalar are likewise similar. The axial and vector examples have an if-factor of around 1.2, which 
corresponds to an increase in A of around 5%. The scale variation for these operators displays 
a similar behaviour to the background process [25] . These processes have no as dependence 
at LO, therefore the only scale dependence enters through the factorization scale dependence 
associated with the PDFs. As a result, the cross section increases with increasing /x which is the 
opposite dependence to the renormalization scale which decreases the cross section as /i increases. 
At NLO one may have expected a large scale dependence for this process, since it is LO in the 
renormalization scale. However, there is a net cancellation between variations in the factorization 
scale and renormalization scales such that the scale dependence for at NLO is very small. The 
axial and vector operators inherit this trait, at LO the (factorization) scale dependence is around 
9% whilst at NLO the scale variation is around 2 — 3%. As is clear from the plots, the variation 
around the default scale at LO does not include the NLO curve so some care should be taken in 
interpreting this small variation as a small theoretical uncertainty. 

It is immediately obvious that the scalar and pseudo-scalar operators do not share the same 
traits as the other operators when it comes to scale variation. This is also naturally explained 
from the discussion in the previous paragraph. Here the heavy sea quarks have a very small 
factorization scale dependence, and as such the LO result barely depends upon the scale choice 
(< 1%). However, at NLO the lack of a factorization scale dependence in a significant part of the 
calculation removes the accidental cancellation. As a result the renormalization scale dependence 
dominates and there is a larger scale dependence at NLO around (10%). As was the case for 
the monojet study the .fT-factor for these operators is larger (1.5-2), however one expects large 
NNLO corrections for these operators since the contributions from two-glue initial states should 
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Figure 7: LO and NLO cross sections for DM production in association with a single photon at the 7 TeV 
LHC. The solid line indicates the cross section obtained with the default scale fi = m xx , the shaded band 
represents the deviation from this scale when the scales are varied by a factor of two in each direction. The 
phase space cuts described in the text (48) have been applied. 



be comparatively large. Note that the A limit for these operators is expected to be very small, due 
to the quark mass suppression, and one expects the full theory to suppress this even further. 

We consider the differential $x spectrum at NLO, as an example we consider the axial operator 
with m x = 100 GeV. Again we compare to the shape of the dominant background (in this case 
Zj) using MCFM to obtain both spectra at parton level. Our results for the spectrum and LO 
to NLO ratios are shown in Fig. [9] and Fig. 10 . Since our explanation of the hardening of the 



spectrum for the signal in the monojet case did not invoke any properties of the recoil object we 
expect the signal to also be harder for this operator. This is indeed what we observe, the spectrum 
for the DM signal is significantly harder than that of the Z7 background. Away from the first bin 
the K- factor is also fairly stable as a function of the $x- The mismatch of the photon px and the 
Px cuts means that at LO there is no contribution to the first bin, since px balance enforces that 
the larger photon px cut is also applied to the $x- However, at NLO there can be events, with 
lower in the first bin. This accounts for the large (actually infinite) K- factor in the first bin. 
At NLO the full fiducial phase space is explored, {fix > 140) resulting in a non-zero cross section 
in the first $x bin. 
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Figure 8: A"-factors for DM production in association with a photon at the 7 TeV LHC. The solid line 
indicates the if -factor obtained with the default scale /i = to x ^, the shaded band represents the Zf-factors 
obtained using scales varied by a factor of two in each direction. The phase space cuts described in the text 
have been applied. 
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Figure 9: NLO predictions for the missing transverse momentum spectrum for signal and background 
obtained using MCFM, for the monophoton process, proceeding via an axial operator. The DM mass is 100 
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Figure 10: The differential if-factor for monophoton production via an axial operator for m x = 100 GeV. 
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Finally we note that care should be taken when estimating scale variations for these processes 
if a jet-veto is applied. Additional information on the topic of scale variation under the application 
of jet-vetos can be found in Refs. [371 SB]- 

5. CONCLUSIONS 

Hadron colliders provide an ideal place to search for DM, offering complementary results to 
those obtained from direct detection experiments. Model independent searches for the pair pro- 
duction of DM require some other visible activity in the event, e.g. jets, photons, or vector bosons. 
We focussed on the cases of monojet and monophoton, for which there are existing searches both at 
the Tevatron and the LHC. Colliders are free from astrophysical uncertainties and are competitive 
with the direct detection experiments for light dark matter and DM with spin-dependent nucleus 
couplings. One of the dominant sources of uncertainty of the experimental results from the LHC 
is systematic uncertainties. This is already true for the monojet searches from the 2011 data set, 
and one assumes that this will continue for 2012 data. In this regard, it is important for the 
theoretical community to aid the experimental effort in any way it can in reducing these uncer- 
tainties. We addressed one such systematic error, the uncertainty associated with the precise rate 
for the signal process. Leading Order predictions suffer from a large dependence on the unphysical 
renormalization and factorization scales. Calculating at NLO accuracy curbs this bad behaviour 
somewhat, with the additional happy byproduct that, since the NLO corrections typically intro- 
duce a iT-factor > 1, limits using the NLO prediction will be stronger in addition to having smaller 
systematic uncertainties. 

To remain as general as possible, and to allow extension of our results to other similar processes, 
we factorized the SM production from the BSM decays. We provided analytic expressions both 
for the SM production at NLO and the LO decay amplitudes of the DM particles. The exact 
nature of the SM production (and DM decay) depends on the particle mediating the interactions, 
and we presented results for interactions which were mediated by vector, axial, scalar or pseudo- 
scalar particles. We considered both the case where the mediating particle was very heavy and the 
effective coupling between DM and the SM was a contact operator and the case of a propagating 
mediator. We implemented all our results into MCFM and they will be made available in the 
next release 2 . This implementation is sufficiently versatile that it can handle all the operators 
considered here, as well as more general (parity violating) variations, as well as more general UV 
completions 

We used our NLO results to study the monojet and monophoton phenomenology in the presence 
of a DM signal. We presented NLO and LO cross sections and demonstrated there is an increase in 
rate as well as a reduction in the renormalization/factorization scale dependence. Thus reducing 
the systematic uncertainty in the signal prediction. Using similar cuts to the experimental analyses 
we were able to infer the amount of improvement on the limits on A obtained by using a NLO 
prediction for the total rate. Typically these limits were increased by around 5 — 20% with the exact 
results being highly operator dependent. Using these NLO results we were able to confirm earlier 
results that saw a significant hardening of the signal $t spectrum, compared to the dominant 
irreducible background. Furthermore, we found a slight enhancement of this feature at NLO due 
to a rising differential i^-factor. 

We also investigated some of the effects of re-instating the mediating particle for the monojet 
case. This confirmed the validity of the EFT for the axial theory. However we saw that the region 
of A probed at the LHC by the top-scalar operator is in a region in which the differences between 



2 Advanced copies can be obtained from the authors, upon request 
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the full theory and the EFT are extremely sensitive to the mediator mass and the exact fiducial 
phase space cuts. As a result we suggest that, should the collaborations proceed to set constraints 
on this operator, they investigate the typical differences between the FT and EFT in order to 
ensure the validity of the EFT results. 

Going forward, as the LHC collects more data at both higher energies and higher instantaneous 
luminosity, it will continue to constrain the dark sector. As running conditions change so must 
analyses and having signal predictions at NLO, through MCFM, will allow the analysis cuts to be 
tuned to best separate signal events from background and maximize the potential for finding DM 
at the LHC. 
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Appendix A: Spinor Helicity Formalism 

In this appendix we define our spinor products, (for a review see ref. [49J) The function u±(ki) 
represents a massless Weyl spinor of momentum k{ and positive or negative chirality. In terms of 
these solutions to the Dirac equation, the spinor products are defined by, 

(ij) = (r\j + ) = u-(k i )u + (k 3 ), (Al) 
[ij] = =u + (k i )u.(k 1 ). (A2) 

We use the convention [ij] = sgn(A;°fc°) (ji)*, so that, 

(ij) [ji] = 2 h ■ kj = . (A3) 

(a\i\b] = (ai) [ib] , (a\(i+j)\b] = (ai)[ib] + (aj)[jb] . (A4) 

Further useful identities are, 

(fVl^) = 2 K , (*l7"tf](*|7/i|*] = 2 (ik) [Ij] . (A5) 

We will also need the following one-loop basis functions, those associated with the reduction of 
tensor triangles, 

t / \ Wx/y) L (x,y) + 1 

L (x,y) = 1 _ x » L l( x >y) = ' ' 

y y 

and the finite part of the one-mass box given by 

LB-if^, ^ =U 2 (l- +U 2 (l- **) + ln (*) In £ (A7) 

\yi J2j V yi/ V J2 J \yij \y2j 6 
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Appendix B: NLO Vector Currents 

We follow the notation of [36J and decompose the amplitude as follows, 



= criAp^V + F») , (Bl) 
where the ubiquitous one-loop prefactor is defined as 

Cr " (4tt)2-« r(l + 2e) ■ (B2) 

The amplitudes we will write down presently have not been UV-renormalised, in order to obtain 
UV finite results one must perform a UV-subtraction, e.g. in the MS-scheme one should subtract, 



cvN cg H-[^-^))A^. (B3) 



2 , If 11 2nf\\ („,„) 
3 3NJ) Av 

We are now in a position to write down the virtual corrections to the monojet process, the leading 
colour contributions have the following form, 

^(1+ 2+ 3Z) = c r (A^ ) V lc + F 1 ^) , (B4) 

with 

^4^y + mi4^y-3, (B5> 



snj \-s 23 J J 2e \ -s 2 3 



and 



(l 9 ,2 g ,3,j- ^ LS -\- Sl23 '- Sl23 ; + 2 (12) (23) 



~L (s 23 /si 23 ) 



1 (13) 2 [1| 7 ^|(2 + 3)|1] U(s 23 /s 123 ) 

+ i — WW) ^~ ■ (B6) 



The subleading in colour amplitude has the following form, 

1 / ^ V 3 / /i 2 V 7 



The remaining _F S ''^ pieces are, 



" 9 ' « y 2 (12) (23) —\- Sl23 '- Sl23/ 
1 (13) 2 (2|(1 + 3)| 7 "|2) /- ai3 -s 23 \ [12] (3|(1 + 2)| 7 ^|1) 

2 (12) 3 (23) V.-S123 -5123/ (12)5123 

1 (1|(2 + 3)| 7 "|1) [12] 2 (23) U(s 123 , s 23 ) , 1 <1|(2 + 3)| 7 "|1) (23) [21] L (s 123 , s 23 ) 



+ 



4 (12) s 2 23 2 (12) 2 S23 

1 (31) [12] (2|7"|2] s l2 3 U(si 23 , s 13 ) , 1 (31) [12] (2|(1 + 3)| 7 "|1> L (si 23 , s 13 ) 



+ 



2 (12) s 2 3 2 (12) 2 sis 

[12][3|(l + 2)| 7 ^|2]-[23] [1| (3 + 2)| 7 "|2] 

4 [13] [23] (12) ' 1 J 
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Finally, we note that for the vector current the diagrams associated with closed fermion loops vanish 
via Furrys theorem, when constructing axial currents we will need the following non- vanishing np 
axial current 

F-»(lt,2+,3T) _ _W|2][21]L l(8l3 ,» m ) 

2 S123 

We note that the above results can be checked by contraction with the current in (4|7^|5] /S45, 
reproducing the formulae listed in |36j . 

At NLO we also require the tree-level amplitudes involving the emission of an additional par- 
ton. The necessary (with the remaining helicity assignments being obtained via line-reversal and 
conjugation) two gluon amplitudes are, 

4 M(1+ 2 + 3+ 4") - m+l±Mm (BIO) 
A v {L q ,2 g ,6 g ,%) - 2 (12) (23) (34) ' [ ' 

A M ( ,+ 2+ o- _ (31) [12] (3|(1 + 2)|7"|4) 
A V [L g ,Z g ,6 g ,^) - 2 (12) s 23 s 123 

(34)[42][2|(3 + 4)| 7 "|1] (3|(1 + 2)| 7 "|(3 + 4)|2] 

2[34] S23S2 34 2(12)[34]s 23 ' 1 ' 

o- 0+ a-, - [13] 2 (2|(1 + 3)| 7 ^|4) 



The four-quark amplitudes are 



2 [12] s 23 s 12 3 
(24) 2 [3|(2 + 4)|7^|1] [13](24)(4| 7 "|1] 
2 (34) s 2 3S 2 34 2 [12] (34) s 23 



+ o+ a- ^ _ (13) (3|(l + 2 + 4)| 7 l3) 
2(12) (14) (23) (34) 



(B12) 



^Ot^W = i(MfcW + MMffl) , (B13) 

H Q ^ q 2 V Sl 2 3«23 S234S23 / 

4°^(l+,2r,3 Q ,4+) = yt^(l+ 23,4+ 3j . (B14) 
For the monophoton calculation we will also need the following tree-level amplitudes, 



45^(1+ 2+ i£,4+) = - X 7^w^wIwo,: ' , (B15) 



4 M (1+ 2+ 3^,47) = ^(1+ 4;,2+ 3r) + ^)(l+ 2+ 4-,3r) . (B16) 
In the above equations we have defined the photon as to emphasize that it is not colour-ordered. 
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